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We discuss three implementations of the CCFM evolution equations in event gen- 
erator programs. We find that some of them are able to describe observables such 
as forward jet rates in DIS at HERA, but only if the so-called consistency con- 
straint is removed. We also find that these results are sensitive to the treatment 
of non-singular terms in the gluon splitting function. 



1 Introduction 

Small- a; final states are central to the understanding of small- a; evolution in 
general. To understand final-state properties such as forward jet rates and 
transverse energy flow, it is important to have event generators which give a 
good description of the experimental data. It has, however, turned out to be 
difficult to produce such an event generator. 

To leading double-log accuracy, the CCFM 0] evolution should be the best 
way of describing small- a; final states, and some attempts to implement CCFM 
into an event generator has been made e.g. SMALLX ^ and LDCMC ||. Below 
we will discuss some recent developments of these generators as well as the 
new CASCADE 0] generator for CCFM. 



2 CCFM and the SMALLX and CASCADE programs 

The implementation of the CCFM Q parton evolution in the forward evolution 
Monte Carlo program SMALLX is described in detail in p| , and we have already 
reported on some recent developments in |0j and . The main new ingredient 
is a modification in the treatment of the non-Sudakov form factor, A„ s , which 
allows for the removal of the so-called consistency constraint. Rather than 
using the simple form, logA ns = — a s \og(l/zi) log(k^ i /ziq^ i ), which requires 
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Figure 1: (a) The cross section for forward jet production as a function of x, compared to 
HI data |7j. (b) The values of z in the initial state cascade for events satisfying the selection 
of forward jet production. The solid (dashed) lines are the prediction of the Monte Carlo 
without (with) applying the "consistency constraint" (c.c). (c) is the same as (a) but using 
the LDCMC program including only gluonic ladders with (dotted line) and without (full line) 
the non-singular terms in the gluon splitting function. 



the consistency constraint > Ziq^ in order to be below unity, the full form [gj 
is used: 

fz \ ( k 2 \ f 1 tf*w/«H>l 

logA„ s = -a s log — log tJ -Y ) ,z = < ku/qu # z < k ti /q ti < 1 (1) 

W ( z if lex/to <z 

This form factor is well behaved for all emissions and gives no suppression in 
the region kti / qu 5: z where then A ns — 1 . 

In Fig. |TJa the prediction for the forward jets is shown. The data are nicely 
described. Including the consistency constraint, the x dependence of the total 
cross section changes, but a similarly good description of F% is obtained by 
changing the infrared cutoff, Qq. However the forward jet cross section is 
reduced. 

The effect of the consistency constraint is well illustrated in the spectrum 
of the values of the splitting variable z for events that satisfy the criteria 
of forward jet production (Fig. [l]b). Only about half of the events satisfy 
the consistency constraint, and especially medium and large values of z are 
rejected. Furthermore we notice that in general z is not very small thus showing 
that even in forward jet production we are far away from the asymptotic region, 
where the small x approximation is valid. 

The CASCADE program which implements CCFM evolution in a backward 
evolution algorithm, gives results which agree well with the ones from SMALLX. 
It should also be noted that these programs also are able to reproduce other 
data, such as F% and photo- production of J /iff as measured at HERA. 



2 



3 The Linked Dipole Chain Model and the LDCMC program 



The Linked Dipole Chain Model g, LDC, is a reformulation of the CCFM 
evolution. The main idea is to reinterpret the non-Sudakov form factor as a 
normal Sudakov for the no-emission probability in some phase-space region. By 
redefining the division between initial- and final-state splittings (which is done 
by angular ordering in CCFM), requiring the transverse momentum of a gluon 
emitted in the initial-state to be larger than the minimum transverse momenta 
of the connecting propagators, qf ± > min(fc^ 1 , k?_ 1± ), an extra weight is given 
to each splitting corresponding to the sum of all emissions which in this way 
are treated as final-state. In the limit of emissions which are strongly ordered 
both in longitudinal and transverse momenta of the propagating gluon, it can 
then be shown that if the kinematical constraint is applied, this sum exactly 
cancels the non-Sudakov. It has, however, not been proven that this holds if 
the consistency constraint is relaxed. 

LDCMC is able to give a good description of F2, but underestimates the 
forward jet rates by approximately a factor 2. Much effort has lately been put 
into the understanding of the differences between SMALLX and CASCADE on 
one hand and LDCMC on the other. Some of this work was reported on in ||. 
One difference between LDCMC and SMALLX is that the latter only includes the 
singular parts of the gluon splitting function, i.e. the l/z and the 1/(1 — z) 
terms, where the l/z term is regularized by the non-Sudakov: 

P g (z) = A ns - + -^—. (2) 

z 1 — z 

LDCMC, however, uses the full splitting function including the non-singular 
terms (where the non-Sudakov has been canceled): 

P g (z) = - + -^—-2 + z(l-z). (3) 

z 1 — z 

In the limits of z close to zero or one, which is the limit in which the CCFM 
equation is derived, the non-singular terms should be negligible. But, as seen 
in fig. ^b, the typical z- values in events with forward jets are around 0.5 where 
the non-singular terms may reduce the splitting function with almost a factor 
2. It is intriguing to note that the forward jet rates obtained by LDCMC is 
approximately a factor two below those obtained with SMALLX (which agrees 
with data). Indeed, preliminary investigations show that removing the non- 
singular terms from LDCMC does increase the forward jet rates as seen in fig. [j]c 
although not with as much as a factor of two M . 
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4 Open questions 

The fact that the results from the SMALLX, CASCADE and LDCMC programs 
depends strongly on the treatment of the consistency constraint and on the 
non-singular terms in the gluon splitting functions, means that also the CCFM 
equation as such is sensitive to these non-leading effects for observables such 
as the forward jet rates. 

To include the non-singular terms in the gluon splitting function in the 
CCFM equation is not completely straight forward. Simply adding them to 
eq. ||may result in negative splitting probabilities. But in the limits of z close 
to zero or one, it is allowed to let the non-Sudakov in eq. || multiply also the 
1/(1 - z) term: P g = A ns P° m 3(z) = A ns (l/z + 1/(1 - z)), in which case 
it is possible to replace the Pg mg (z) with the full splitting function of eq. |[ 
Preliminary investigations show that doing this in the SMALLX program 
indeed gives a large effect, although more studies are needed to quantify the 
influence on the results for e.g. forward jet rates. 

In conclusion, it is still an open question whether the CCFM evolution 
equation is an appropriate way of describing small- x final states or not. Much 
more work is needed, both phenomenological work using the event generators 
and purely theoretical studies of the CCFM evolution, to understand how 
to correctly treat the non-leading effects which seem to be important in the 
description of data. 
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